Nonlinear model for magnetosonic shocklets in plasmas 
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Abstract 

Exact nonlinear equations for magnetosonic shocklets in a uniform hot magnetoplasma are de- 
rived by using the nonlinear magnetohydrodynamic equations. Analytic as well as numerical solu- 
tions of the nonlinear equations are presented. Shock-like structures of the ion fluid velocity and 
magnetic field (or the plasma density) perturbations are obtained. The results may have relevance 
to the understanding of fast magnetosonic shocklets that have been recently observed by onboard 
instruments of the Cluster spacecraft at the Earth's bow shock. 
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In the past, there has been a great deal of theoretical interest (e.g. Refs. |l[ Q, Q]) in 
studying the properties of nonlinear magnetohydrodynamic waves in plasmas. It has been 
found that both slow and fast magnetosonic (FMS) waves can propagate in the form of 



either solitary or shock waves in plasmas. Very recently Stasiewicz et al. jj, ^ reported 
detailed properties of slow magnetosonic (SM) solitons 4| and FMS shocklets P|, which 
have been observed by a fleet of four Cluster spacecraft at the quasi-parallel bow shock. 



Observations reveal that SM solitons 



a 



are associated with large amplitude compressional 
(rarefactional) magnetic field (plasma density) variations. On the other hand, FMS shocklets 
are accompanied with compression of the plasma density and magnetic field perturbations. 

In this Brief Communication, we present a nonlinear model for FMS shocklets, which 
may account for the observed FMS shocklets at the Earth's bow shock. Specifically, we 
show that FMS shocklets are associated with the nonlinear steepening j?] of arbitrary large 
amplitude FMS waves in a high-beta magnetoplasma. 

The dynamics of the nonlinear FMS waves in a magnetized plasma is governed by the 
inertialess electron momentum equation 

= —rieeE — S/pe — Uee— x B, (1) 

c 

where rig is the electron density, e is the magnitude of the electron charge, E is the wave 
electric field, B is the sum of the ambient and wave magnetic fields, Pe = UeTf, is the electron 
pressure, Tg is the electron temperature, Ve is the electron fluid velocity, and c is the speed 
of light in vacuum. The ion dynamics is governed by the ion continuity equation 

^ + V ■ (n,v,) = 0, (2) 

and the ion momentum equation 

p + Vj ■ V j Vj = UieE - Vpi + Uie^ x B, (3) 

where is the ion number density, Vj is the ion fluid velocity, p = niirii is the ion mass 
density, mi is the ion mass, pi = riiTi is the ion pressure, and Tj is the ion temperature. 
Equations (l)-(3) are closed by means of Ampere's and Faraday's laws 

Attc 

V X B = {riiVi - neVe), (4) 
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^ = -cVxE, (5) 

together with the quasi-neutrahty condition Ue = rii = n. The latter is vahd for a dense 
plasma in which the ion plasma frequency is much larger than the ion gyrofrequency. Equa- 
tion (4) holds for the FMS waves whose phase speed is much smaller than the speed of 
hght. 

Eliminating E from (1) and (3) we obtain 

dvi „ (V X B) X B ^2^, 

where the ion sound speed is denoted by Cg — [{T^ + Ti)/miY^'^. On the other hand, from 
(1), (4) and (5) we have 



dB _ 



Vi X B ^(V X B) X B 



(7) 



We are interested in studying the nonlinear properties of one-dimensional FMS waves 
across the external magnetic field direction zSq, where z is the unit vector along the z 
axis and Bq is the strength of the ambient magnetic field. Thus, we have V = x9/9x, 
Vj = MX and B = B{x)z, where x is the unit vector along the x axis in the Cartesian 
coordinate. Normalizing p by the equilibrium mass density po — nouii, u by the Alfven speed 
Va — Bo/y/ATTpo, B by Bq, time by the ion gyroperiod x by Va/oJci — c/uj^i, where uj^i 
is the ion plasma frequency, we have our nonlinear MHD equations in dimensionless form as 

and 

dh ^du dh ^ 

dt dx dx ' 
where h = B{x,t)/ Bq and f3 = C^/V^ = 47rno(Te + Tij/B^ represents the plasma beta. 

Equations (8) and (10) yield p = h, a. concept of frozen- in- field lines in a magnetized 
plasma. Hence, we have from (9) and (10) 
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and 

dh ^ dh ^ ^du Q 
dt dx dx 

In the zero-/5 limit, Eqs. (11) and (12) agree completely with Eqs. (2a) and (2b) of Stenflo et 
al. jsl who demonstrated rapid steepening of the velocity and magnetic field perturbations 
leading to the formation of FMS shocklets in a cold magnetoplasma. 

In the following, we study the properties of FMS shocklets in a warm magnetoplasma. 
We introduce the change of variables 



..2V^.>„(^^).c.. (13, 



where ci and C2 are constants. This diagonalizes the system of equations (|TTj) - (|T2j) to the 
form 



^'*' + («. + V^)?l = 0, (15) 



dt dx 

where u and h are given in terms of ipi and il)2 by Eqs. (fT!?|l - (fT^ . This system admits 
"simple wave" solutions [9j, which can be found by either setting ipi or ip2 to zero. Setting 
ip2 to zero in Eqs. (fT3 j) - (fTB|) . we obtain 

t^e"+y^^)t^o, (17) 

where from Eqs. (|T!^ - ()14|1 



4 



The constant 



is determined by the boundary conditions u = and h = 1 at |x| = oo. Using ipi = 2u in 
Eq. (fTTjl . we have 

|^(.,y^)g.o, (19) 

where h{u) is imphcitly given by Eq. ()18|) . We note that the results given by Eqs. fllSp 
and (|T9|) generahze the resuhs presented in Ref. for the case of arbitrary /5 values. The 
paths in the {x,t) space where u is constant, can be described by the ordinary differential 
equation 

f = «+V^' (20) 

where the right-hand side is constant (since u is constant along the path), which after 
integration gives 

x{t) = {u + ^Jh + (3)t + Xo. (21) 

Here xq is a constant of integration. The general solution of Eq. (fT^ is a function of the 
integration constant xq, viz. 

u = /(xo) = fix - (m + ,/h^)t], (22) 

where h{u) is given by Eq. (|18|) and / is a function of one variable, determined from the 
initial condition at t = 0; the velocity u can be evaluated for different x and t by solving 
Eq. for u. Equation (j^^ describes a nonlinear FMS wave propagating in the positive x 
direction, where the time-dependent solution has a typical structure of the wave-steepening, 
similarly to the solution of the inviscid Burgers equation This solution may also be 
obtained by directly assuming that h can be written as a function of u. Note that there are 
no steady state solutions within this system, since the dispersion is absent; but by including 
the electron inertial effect in Eq. (fTUI) this could be achieved on length scales < Ag = c/upe, 
where Upe is the electron plasma frequency. The dispersive effects break the p = h relation, 
and produce an asymmetry between the density and magnetic field perturbations. 
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We have analyzed the system ()lH ) - p2p numerically. As an initial condition, we took the 
magnetic field h = 1 + 0.5 sech(a;/20), describing a localized magnetic field (and density) 
compression of the plasma. As we are interested in the evolution and creation of a "shocklet" 
moving in one direction, we choose the initial condition for the velocity u from Eq. (|18|). 
describing a wave moving in the positive x-direction only. The evolution of the system for a 
low-beta (/? = 0) and for a high-beta (/? = 10) plasma is displayed in Figs. 1 and 2, respec- 
tively. In agreement with the analytical prediction, the initial wave is propagating in the 
positive x-direction only. Figure 1 shows that both the magnetic and velocity perturbations 
steepen and a shock front starts to develop. For the high-beta plasma displayed in 

Fig. 2, the velocity associated with the shocklet is larger, and the self-steepening develops 
faster than for the low-beta plasma case. When the shock fronts become steep enough, ef- 
fects such as the electron inertia and electron Landau damping (which heats the plasma) will 
become important. Dispersive effects also occur if the waves propagate with some angle to 
the magnetic field direction. The combined effects of dispersion and wave-particle induced 
dissipation could explain the apparent phase asymmetry between the magnetic field and 
density perturbations, as observed in large-amplitude FMS shocklets This asymmetry 
is likely to appear after the shocklets have developed due to the self-steepening of the FMS 
waves, as investigated here. 

To summarize, we have considered the nonlinear propagation of FMS waves in a hot 
magnetoplasma. It has been shown that the nonlinear MHD equations in a finite-/? plasma 
can be reduced to a pair of equations in which the ion fiuid velocity and the compressional 
magnetic field are nonlinearly coupled. The system has been diagonalized and special, single 
wave solutions have been obtained. The solutions represent the spatio-temporal evolution of 
an arbitrary amplitude FMS waves. The equations for the full system are solved numerically 
to show the formation of FMS shocklets, in full agreement with the analytic results. The 
finite plasma beta has significant infiuence on the shocklet profile in that the shocks develop 
on a much shorter timescale than for the the plasma with a low beta. In conclusion, the 
present results qualitatively account for the salient features of the observed FMS shocklets 
at the Earth's bow shock Sf. 
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FIGURE CAPTIONS 



FIG. 1. The evolution of the normahzed ion fluid velocity (upper panel) and compressional 
magnetic field (lower panel) for /3 — 0, for the times 1) t — 0, 2) t — 12.25 3) t — 24.75 and 
4) t = 37.25. 

FIG. 2. The evolution of the normalized ion fluid velocity (upper panel) and compressional 
magnetic fleld (lower panel) for P — 10, for the times 1) t — 0, 2) t — 12.25 and 3) t — 24.75. 
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